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An attempt is made to assess the significance of rotation in the core- 
collapse supernova phenomenon, from both observational and theoretical 
point of view. The data on supernovae particularly indicative of the role 
of rotation in the collapse-triggered explosion is emphasized. The problem 
of including the rotation of presupernova core into the supernova theory 
is considered. A two-dimensional classification scheme of core-collapse su- 
pernovae is proposed which unifies "classical" supernovae of type Ib/c and 
type II, "hypernovae" and some GRB events. 

PACS numbers: 97.10.Kc, 97.60.-s, 97.60.Bw, 97.60.Jd 



1. Introduction 

The phenomenon of supernova amazed already ancient observes as some 
bright historical supernovae were visible on the sky even in the daytime. The 
absolute luminosity of supernovae has been properly estimated only in the 
20-th century, when it was realized that supernovae belong to a very special 
class of astronomical events. In 1885 the nova S And appeared in the M31 
nebula, now well-known as the Large Galaxy in Andromeda. At those times 
many astronomers accepted the in-Galaxy theory of M31 and other nebulae. 
After establishing that the real location of M31 is extragalactic, astronomers 
were forced to conclude the nova S And^ was much brighter than any usual 
nova ^ - it was a super-nova! 



^ odrzywolek@th.if.uj.edu.pl 

^ Now called SN 1885A. 
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The systematic supernova research began in the 20-th century. Unfor- 
tunately, there was no Galactic supernova event since the 17-th century. In 
spite of this astronomers have observed more than 2000 extragalactic super- 
novae. The number of observed events grows rapidly, from about 20 per year 
in the eighties to about 200 per year now. In contrast to the optical events, 
more than 600 supernova remnants have been found in the Galaxy. Also, 
a number of extragalactic remnants, mainly in the Local Group galaxies 
LMC, SMC, M31 and M33 have been observed. 

In 1942 Minkowski |3j introduced the modern classification scheme of su- 
pernova events into two classes. To the first class belong supernovae with no 
hydrogen absorption lines in the spectrum referred to as type I supernovae. 
The second class comprises the supernova events with strong hydrogen lines 
which are referred to as type II supernovae. 

As for the physical nature of supernovae. Landau Q in 1932, soon after 
discovery of the neutron, suggested the possibility of existence of dense 
stars composed of neutrons which are stabilized by very high pressure of 
the neutron gas. In 1939 Baade&Zwicky [5j proposed the gravitational 
collapse of a normal star to such a neutron star as the supernova energy 
source. This picture is generally accepted today. In 1960 Fowler and Hoyle 
[HI [Tj pointed out that also nuclear reactions can serve as a source of the 
supernova energy. They proposed the thermonuclear explosion of a white 
dwarf or a giant star as an another possible supernova mechanism. 

After several decades that passed since the original proposal of Zwicky 
and Baade, significant progress in understanding the supernova mechanism 
has been achieved. The detection of the neutrino burst correlated with the 
appearance of SN1987A proved that the theoretical research is on a right 
track. Unfortunately, the standard core-collapse supernova theory suffers 
from difficulties in producing a successful explosion under general condi- 
tions. This could be a consequence of suppressing star's rotation in the 
theory. In this paper we collect arguments in favour of the necessity to in- 
clude rotation into the supernova theory. The paper is organized as follows. 
In the next section we review in some details the modern version of super- 
nova classification. We wish to emphasize a division of observational data 
into these related to outer layers of the exploding star and the ones reflect- 
ing the physics of the engine mechanism. In Sect 15] the essential features 
of the standard supernova theory are reminded. The problem of inclusion 
of the presupernova core rotation into the supernova theory is considered 
in SectJlJ Finally, in Sect El we introduce a two-dimensional classification 
employing some measure of the rotation as a second dimension. 
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2. Modern classification of supernovae 

With a simple classification of supernova events available, researchers 
were tempted to establish a direct relationship between the two proposed su- 
pernova mechanisms and the two observed supernova types. Unfortunately, 
the relation of the thermonuclear explosion mechanism and the core-collapse 
scenario to type I and type II supernovae is not as straightforward as it may 
seem. It turns out that type I supernovae with no hydrogen lines form a 
class of diverse events which can be further divided into more homogeneous 
subclasses. The spectra near the maximum brightness are used to distin- 
guish type la events which show a strong Sill absorption dip at A 6150 A 
and type Ib/c events with no or weak Sill lines. Further, a strong helium 
line at A 5876 A is employed to distinguish type lb events from those of 
type Ic. Type Ib/c supernovae occur in the same environment as do type 
II supernovae, namely in the star forming regions in spiral galaxies. Both 
are thus related to young population I stars. Type la supernovae are more 
common and they occur in all types of elliptical and spiral galaxies, and in 
the halo of our Milky Way galaxy. This location indicates they are related 
to population II stars. The main observational features used as a basis of 
modern supernova classification are presented in Table I . 

A generally accepted hypothesis with regard to the relationship of the 
physical supernova mechanism and the observed supernova types postulates 
that thermonuclear explosions of white dwarfs produce type la supernovae 
and the core collapse of giant stars results in type Ib/c/II supernovae. Basic 
arguments in favour of this identification come from common locations of 
type la supernovae in old stellar systems and type Ib/c/II - in regions of 
active stellar formation in spiral galaxies. 

Homogeneity of type la supernovae also indicates the same explosion 
mechanism. Type la events have similar spectral properties and obey a 
simple empirical relation {e.g. the Phillips relation 0) between the maxi- 
mum absolute brightness and the behavior of the light curve, which states 
that brighter events are longer. The observed scattering of absolute mag- 
nitudes at the maximum is about 2™, and further division of type la into 
subclasses is not excluded [H]. 

It was long believed that these differences are due to observational errors 
and all type la supernovae are almost identical. However, the maximum 
brightness and the decay time are possibly related to the amount of ^^Ni 
synthesized during the explosion. The light curves clearly show that the 
beta decay of ^^Ni into ^®Co, and a subsequent decay of ^^Co to ^®Fe with 
a half-life of 77 days, are the main energy sources supporting the light 
emission during the late phase of the explosion. It is of great interest that 
the radioactive energy source from ^^Ni and ^^Co decays occurs in all types 
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Table I 



Supernovae classification and properties. 



Type 


la 




II 


Definition 


No hydrogen lines 


Strong 
hydrogen lines 


Location 


All type galaxies and 
halo 


Spiral arms, star forming regions 


Relative rate 

Galaxy rate 


30% ± 9% 
~1/200yr 


11%±6% 
~1/500yr 


59% ± 28% 
~1/100yr 


Absolute 
magnitude 


-20™ 
homogeneous" 


(-19.5™)* 


-16™. -19™ 
(-lO"")* 


Spectra at 
maximum 


Si II absorption 
at X6355A 


No or very 
weak Si W' 


Strong Ha 


Late spectra 

(months 
after max.) 


[Fen]+[Fein] emission 


Forbidden [OI] + [Call] emission^ 


Nebular 
remains 


Shell 
(B aimer-dominated) 


Shell 
(Oxygen-rich) 


Plerion 
(Crab-like) 


Compact 
remains 


None 


Neutron star ( or a black hole ) 


Progenitor 


Accreting white dwarf in 
binary system 


Wolf-Rayet star 


Giant star 


Main light 

curve energy 
source 


^^Ni 6 days ^^^Co -> 77 days ^^^Fe decay 


Explosion 
enei^y source 


Fusion C/O — > Fe 


Neutrino flux 


Physical 
mechanism 


White dwarf 
thermonuclear disruption 


Core-collapse 



"Average Mb ~ -19.8±0.03 + 51og(Ho/75), but deviations do exist: SN1991bg -le.S", SN1991T 
-20.5"". Ho is Hubble const. [km/sMpc ']- 
*Extremly bright events called hypemovae. 

"^Doublet >.6347A,^637lA blueshifted (velocity up to 30000 km s"') to ~X6150A. 
''Absorption of Hel (X5876A) defines type lb; lack of Hel - type Ic. 
"Ha always dominates in type n SNe. 
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of supernovae in spite of very different explosion mechanisms. 

Significant constraints on the supernova theory are imposed by observa- 
tions of the nebular remnants. Unfortunately, in spite of known positions 
of SN1006, SN1054, SN1572 and SN1604 we can not firmly classify these 
historical supernovae. In effect, it is not perfectly clear what remains after 
type la/b/c and II supernova explosions. We can distinguish three basic 
types of remnants: 

• Balmer-dominated shell 

• oxygen-rich shell 

• plerion (Crab-like nebula) 

A classical example of a plerion is Ml, the Crab nebula. Because the 
pulsar is found in its center we connect plerions to type II events. They 
are sometimes surrounded by weak shells. The Balmer-dominated shell is 
usually connected to type la, and its emission can be explained as a result 
of the interaction between the shock wave and the interstellar hydrogen. 
Arguments confirming this identification come from the analysis of SN1572 
and SN1604 remnants, which both probably were type la supernovae. In 
addition, we can not see any compact stellar remnant, such as a pulsar or a 
neutron star, inside or near outside the shell. Finally, type Ib/c supernovae 
are supposed to leave as a nebular remnant the oxygen-rich shell. This 
identification is supported by the theory of origin of type Ib/c supernovae 
from stars stripped off hydrogen and helium layers. The lost matter at 
least partially fills in the space around donor stars. If extensive mixing 
and core dredge-up took place at some evolution phases of these stars, we 
expect some enrichment of circumstellar medium (CSM) in heavy elements. 
The shock interacts with CSM whose composition is a relic of star's earlier 
evolution and excites atoms of heavier elements, like e.g. oxygen. 

Observed spectral differences between various types of core-collapse su- 
pernovae are in significant part a result of the differences in the outer shells. 
Classification of core-collapse events clarifies and becomes more elegant 
when we use the lost mass as a basic ordering parameter |1U| . Stars with 
very massive hydrogen layers at the onset of the explosion produce type 
II-P supernovae with the light curve plateau. An extreme example of such 
a light curve was provided by the SN1987A. If the star has not very massive 
hydrogen shell the explosion is classified as type II-L with fast, linear decay 
of the light curve. In case of a very thin H shell the explosion is type lib and 
if the whole hydrogen is lost, the event is classified as lb. If in addition the 
helium layer is also removed we can see type Ic supernova. This ordering of 
core-collapse supernovae is displayed schematically in Fig. ^ 
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II-P II-L lib lb Ic 

SN1987A SN1993J SN1998bw 

Very low Medium High Extreme Mass loss rate 

5-10 Mq 1-2 Mq -0.2 M0 NoH, No H Outer 

H shell H shell H shell massive He & no He shells 

i — stellar wind, binary mass exchange? — > 



Fig. 1. Unification of the core-collapse events within an one-dimensional family of 
the growing mass loss rate. 



Within this scheme, every type of core-collapse supernovae covers a 
range of outer shell masses of pre-supernova stars. For every pair of su- 
pernovae we may find a third one which could be placed between them 
in the classification scheme. For example, SN1993J in M81 was a missing 
link between type II and type lb with a residual hydrogen shell, now distin- 
guished as type lib. The evidence of supernovae corresponding to stars with 
thin He layers is not clear yet, but candidates exist The fact this clas- 
sification works indicates that mass loss rates in progenitors increase from 
type II to type Ic. Outer layers are expelled by strong stellar winds, but 
in extreme cases also stripping due to the mass exchange in close binaries 
occurs, as we indicate in Fig. ^ 

This scheme clearly shows that the engine mechanism is not of major 
importance for the classification of supernovae based on light curves and 
spectra, which both are formed by processes occurring in outer layers of 
the exploding star. There are, however, a number of observations related 
to supernovae, such as pulsar initial velocities, asymmetry of explosions, 
"hyp er novae" , jets and possibly some gamma-ray bursts (GRB), which are 
likely related to processes in the inner part of the exploding star, in partic- 
ular they can directly reflect physical conditions prevailing in the collapsing 
core. The detailed picture of the collapse depends on the mass of the core 
related to the ZAMS mass and the evolutionary track of the progenitor. In 
this paper we focus on the rotation of the progenitor core which can very 
significantly affect physics of the core collapse. Below, in sectE) we discuss 
interpretations of the above observations as an evidence of rotation of the 
collapsing core. We propose a two-dimensional classification of the core- 
collapse supernovae which can accommodate "hypernovae" and GRB's. It 
is an extension of the scheme shown in Fig. Q with the second dimension 
being some measure of the amount of rotation of the collapsing core. Fig. |1] 
on page 
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Table II 

The evolution of 25Mq, Z = 0.02 star according to ^21- Duration of a given 
phase is defined here as time since central ignition of a given reaction until 
the central ignition of a next reaction. This includes phases of off-center 
burning, with no reactions in the center. The temperature and density 
correspond to the beginning of each phase. 



Burning 


T,[K] 


Pc [g/cm^] 


Duration 


H 


3.7 X 10^ 


3.8 


5.8 mln yrs 


He 


1.6 X 10** 


200 


85 000 yrs 


C 


6.3 X 10** 


10^ 


280 yrs 


Ne 


1.2 X lO** 


2x 10"^ 


300 days 


O 


1.7 X lO** 


4x lO*" 


134 days 


Si 


2.8 X 10** 


3.2 X 10^ 


30 hours 


Shell Si 


3.4 X 10** 


3.2 X 10^ 


5.5 hours 


Core 
collapse 


1.6 X 10** 


1.6 X lO'' 


0.1 ...0.5 sec 



3. Standard Supernova Mechanism 

In the standard scenario of the core-collapse supernova without rota- 
tion, a massive star at the onset of collapse is a red or blue super-giant, for 
example, the progenitor of SN1987A was a B3I spectral type star. It con- 
tains a number of shells, which correspond to the following stages of nuclear 
burning, given here in a simplified form: 
H — > '^He 

4jjg y 12q 16q 

( ) 20^;, 24Mg 

20Ne ( 160, 24Mg ) 160, 24Mg 

160 ( 24Mg, 28Si ) 28si 
28si(32g) _^ 56^1, 56pg^ 54pg 

The synthesis of the most strongly bound nuclei near ^^Fe ends the 
network of thermonuclear reactions. In Table II we show main evolutionary 
phases of a star of 25Mq. One can notice that burning of neon begins about 
one year before the star's death. 

The late phases of nuclear burning are quite complicated and still not 
perfectly understood. Soon after ignition of an off-central silicon burning the 
core composed of iron-group nuclei loses its stability and starts to collapse. 
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A short explanation is that the core mass exceeds the Chandrasekhar Hmit. 
The full expression of the Chandrasekhar mass is |13j : 

Mch = 1.44 Mo (2 Fe)' 

(1) 

It is a function of the electron fraction Yg) the entropy per baryon 5e, the 
nuclear composition expressed through an average nuclear charge Z, and 
the ratio of radiation and matter pressure Prad/Pmat] ol ~ 1/137 is the fine 
structure constant. For Y^. = 0.5, 5e = and Prad = we obtain the famous 
result = 1.44 Mq. The mass of the stellar "Fe" core is in the range 

1.2 - 2 Mq. a number of processes are important in the last day before 
the collapse. The electron capture by iron and nickel isotopes decreases the 
electron fraction in processes such as e.g. : 

^^Fe + e" — Mn + i^e. 

Neutrinos carry away the energy from the core decreasing its entropy. The 
silicon burning Si — > Fe stops in the core and begins in a surrounding shell. 
The most important reason of the collapse onset is the decrease of entropy. 
Usually the Chandrasekhar mass limit is not really exceeded because of a 
huge external pressure not included in Mch- The instant when the core 
starts to collapse is deduced from the evolutionary code results. When the 
speed of contracting matter exceeds some in advance prescribed value, the 
evolutionary track is finished. This moment is believed to be close to the 
real stability loss and to start of actual collapse of the core. 

Once started, the collapse is fast, close to a free fall. In the numerical 
simulations the collapse lasts 100 - 500 milliseconds, depending on the initial 
conditions and on the input physics. During the first phase the collapse is 
homologous - the velocity is proportional to the radius. The sound speed 
decreases with the radius, and at a sufficiently high distance from the center 
it falls below the matter speed. This place is called the sonic point. Time 
when the central density reaches its maximum defines the core bounce. At 
this time the infalling matter is stopped by the pressure component due 
to nuclear forces which starts to grow very rapidly. In milliseconds the 
velocity goes to zero and simultaneously the density grows by a few orders 
of magnitude ^[]. Strong sound waves start to propagate outside. As 
a result of non-linear hydrodynamical effects and the supersonic speed of 
matter discontinuity forms in velocity, density, pressure and entropy ~ the 
sound wave becomes the shock wave. 

It was accepted for a long time that the sketched above process is a main 
part of the real explosion mechanism and is able to provide at least a few 



i + fA 
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foe^ of energy. The shock was beheved to traverse the entire star and reach 
the surface to produce a spectacular supernova type II (or Ib/c) event. This 
scenario is known as a prompt explosion. Unfortunately, more detailed cal- 
culations with improved physical ingredients, such as a realistic equation of 
state (EOS), general relativity corrections, sophisticated progenitor models 
and a complete set of nuclear reactions included in new gas-dynamical al- 
gorithms, have shown the failure of this idea. The shock wave gets stuck in 
the Fe core as a result of the energy loss chiefly due to heavy nuclei disso- 
ciation (~ 8.7MeV/nucleon ~ 1.7foe/O.lM0) and - when shock moves out 
of the neutrino trapping sphere - the neutrino emission. Although it is still 
possible to produce a successful prompt shock under extreme assumptions 
on the equation of state or other part of included physics, nowadays neu- 
trino processes are believed to play a major role in core-collapse supernova 
explosions, so let's look at them. 

As we pointed out already, beyond the nuclear matter density, the EOS 
stiffens very rapidly. Simultaneously matter initially composed of nuclei 
and electrons, transforms into a nucleon-electron gas and finally becomes 
an almost pure neutron matter. Other phases of dense matter, such as e.g. 
the kaon-condensed nuclear matter or the quark-gluon plasma can possibly 
form, but usually are not considered in the standard supernova scenario. Al- 
most every electron-proton pair is transformed into neutron and neutrino. 
Initially neutrinos escape freely. But as a result of neutrino cross-sections 
growth with temperature and density, neutrinos become trapped. It means 
that the diffusion time is bigger than the dynamical time-scale. The grav- 
itational energy released in the collapse is "frozen" as the energy of the 
lepton Fermi sea. The inner part of the star which is not blown off, after 
the core bounce and shock traversal becomes a protoneutron star (PNS). 
The edge of the neutrino-trapped zone is referred to as the neutrinosphere, 
with analogy to a photosphere of normal stars. Definition of the neutri- 
nosphere is somewhat ambiguous as different neutrino flavors have different 
and diffused neutrinospheres. Harder neutrinos escape from bigger radius 
than the do soft ones, because of the cross-section dependence on the neu- 
trino energy. The most important fact is the location of the neutrinosphere 
between the shock and the neutrino-rich protoneutron star. The trapped 
neutrinos diffuse out of the protoneutron star in time of the order of one 
second, and simultaneously the protoneutron star contracts from an initial 
radius of about 60km to final radius of 20km, like a Kelvin-Helmholtz star 
jl5j . Neutrinos carry away about 100 foe of energy, i.e. ~99% of the total 
energy released in the collapse. The neutrino flux vanishes in tens of sec- 
onds. This picture was in general confirmed by the detection of neutrinos 



^ foe - fifty-one-erg - lO^'^ erg 
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during the SN1987A event and now this is the most well established part of 
core-collapse supernova theory. In spite of very small neutrino cross-sections 
about 1% of the energy is transferred to the hot radiation bubble between 
the nascent neutron star and the shock in hundred milliseconds. It causes 
the shock wave revival. Later, the explosion is similar to that of a prompt 
mechanism. This scenario is called a delayed or neutrino-driven supernova 
mechanism. Currently it is a subject of very active studies. The neutrino 
transport and convection appear to be the most important processes of this 
scenario [THj . 

The shock produced in the engine of supernova by prompt or delayed 
mechanism traverses the entire star ionizing the gas, igniting nuclear re- 
actions and triggering convection jl7l I18j . The shock wave may reflect at 
the shell boundaries. The onion-like structure of different nuclei layers is 
destroyed. Intensively mixed matter from the center may move close to the 
surface. In a few hours the shock reaches the surface. The photosphere 
begins to expand and we can see the enormous growth of brightness. A 
supernova appears on the sky. During the first few weeks the light curve 
is dominated by the recombination of ionized atoms. Later, the light curve 
mimics the decay curve of beta-radioactive ^^Ni. The subsequent decay of 
^^Co with Ti/2 = 77 days is the energy source for supernova during next 
months. 

4. Rotation in Supernova Theory 

A major deficiency of the sketched above supernova theory is suppres- 
sion of rotation of the presupernova star. Prom the observational point of 
view there is growing evidence that rotation can play an important role in 
the explosion. Below we discuss some relevant data. However, before we 
attempt to asses if inclusion of rotation is necessary to explain the observed 
properties, we address the following important questions: 

• are massive stars rapid rotators? 

• do cores of massive star rotate? 

• what is the shape of the rotating core? 

• what is caused by the angular momentum conservation during col- 
lapse? 

4.1- The simplest rotating star 

Let us consider the behavior of a rigidly rotating body which is incom- 
pressible and homogeneous and bound by the Newtonian gravity. This is 



Core- Collapse Supernova Mechanism - Importance of Rotation 



11 



very simple and idealized model of a rotating star, closer to the liquid drop 
model of the atomic nucleus except of the surface tension, which is 
negligible for astrophysical objects and huge for atomic nuclei.^ It illus- 
trates, however, basic features of the rotation's influence on properties of 
real objects. 

The most important parameter, which determines physics of rotating 
and gravitating bodies, is the ratio of rotational to gravitational energy, 
Erot/Egrav (here Egrav is the absolute value of the gravitational energy). 
In case of zero angular momentum the body is of spherical shape. Any 
rotation leads to a spheroidal shape, as shown by Mac Lauren in 1742. But 
when Erot/ Egrav exceeds the value of 0.1375 two solutions of the problem 
exist (Jacobi, 1834). The first solution is the Mac Lauren spheroid, and 
the second one is a triaxial ellipsoid. The latter is the ground state corre- 
sponding to the minimum of the sum of rotational and gravitational energy. 
Transition from the Mac Lauren spheroid to the Jacobi ellipsoid requires 
some dissipation mechanism, because of difference in the total macroscopic 
energy. Mac Lauren spheroids beyond Ej-ot/ Egrav = 0.1375 are secularly 
unstable with respect to dissipative processes, such as, e.g. viscosity. 

When the rotational frequency of the body increases, for example in 
effect of shrinking, two cases are possible. If the dissipative time-scale is 
short compared to the dynamical time-scale, the body evolves through a 
sequence of Jacobi ellipsoids. For Erot/ Egrav > 0.16 they are unstable secu- 
larly and form so-called Poincare pears. A bit more rotational energy leads 
to the dynamical instability and to the fission. If dissipative processes are 
slow the body evolves along a sequence of Mac Lauren spheroids, which are 
dynamically stable up to Erot /Egrav = 0.2738. Beyond this limit no stable 
configuration exists and the body has to get rid of the angular momentum 
or to breakup. 

Surprisingly, inclusion of a differential rotation and a compressible equa- 
tion of state (and also general relativity corrections) results in minor changes 
of the secular and dynamical instability limits of Erot / Egrav which are im- 
portant for supernova theory. Recent publications report the possibility of 
significant decrease of those limits, for a toroidal density stratification. The 
secular instability limit of the ratio Erot /Egrav becomes 0.038 [201; tran- 
sitions leading to bar- like configurations. The dynamical instability limit is 

found to be Erot/ Egrav = 0.14 |21) . 

During the collapse of the stellar core the ratio Erot / Egrav grows signif- 
icantly, by a factor of a few tens. Given the insensitivity of the instability 
limits with respect to details of rotation and to the compressibility of mat- 
ter, the values listed above seem to be sufficiently accurate to allow us to 

In a ball of water the gravitational energy exceeds the surface tension energy if radius 
is bigger than 10 m. 



12 



A. Odrzywolek et al. 



assess the possibility of triaxial deformations and the core breakup in SN 
events. 

Rotation of supernovae progenitors 

Are supernova progenitors rapid rotators? The answer is well known to 
astronomers Single stars which are supposed to be supernova progen- 
itors begin their lives as O and B main sequence stars with initial masses 
M > 8Mq. The surface velocity is very high for these stars. It can be de- 
termined observationally with some uncertainty due to an unknown angle 
between the rotation axis and the observer direction. The velocity may be 
close (~70%) to keplerian velocity at the surface radius. But, except of the 
Sun^, it is impossible to determine rotation inside the star. To address the 
problem one must resort to stellar modeling. 

In last few years new detailed evolutionary calculations of rotating stars 
have been carried out pil. The results are in good agreement with ob- 
served surface properties. One can thus treat predictions with regard to 
the internal structure with some confidence. Contrary to previous opin- 
ions, numerical results show that cores of stars rotate fast with velocity not 
significantly dependent on the initial conditions. Nevertheless, one should 
mention that the calculations described above neglect magnetic fields which 
may transport angular momentum and slow down the core. 

Typically, according to 521] Fe cores have radii of about 2000 km, masses 
of LSA/q and rotate differentially with periods of the order of 10 sec- 
onds. Outer convective shells surrounding cores rotate rigidly, and the 
rotational frequency drops at shell boundaries discontinuously. The ra- 
tio Erot/ Egrav reaches near the center the maximum values up to 0.04. It 
was assumed in calculations of fSlj that it is safe to neglect the possibility 
of non-axisymmetric core deformations because Erot/Egrav = 0.04 is only 
~30% of the Mac Laurin spheroid secular instability limit. Results cited in 
the previous subsection 20 suggest the secular instability for toroidal "fiz- 
zler"^ configurations at 0.038 and the possibility of "triaxial" cores during 
last phases of the stellar evolution seems not to be excluded if the core is 
far from a spheroidal shape. In spite of the rapid core rotation almost all 
angular momentum is located in distant massive shells, but Erot/Egrav is 
very small there, as indicated in Table III. 

* Helioseismology allows us to see the interior of the Sun. Other stars are in the range 

of observational abilities too; see e.g. . |28) . 
^ This is a transient object between those stabilized by the degenerate electron gas 

pressure {e.g. "Fe" cores and white dwarfs) and neutron stars. Its existence and 

dynamical stability are results of extremely strong centrifugal force, but "fizzlers" 

are secularly unstable. 
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Table III 

Parameters of rotation of inner shells of a typical presupernova star 
from 



Part of 
star 


Angular 
velocity 


Erot j Egrav 


Angular 
momentum 


"Fe" 
core 


10. ..0.01 

rad/s 


0.04... 0.01 


2 X lO'^^erg x s 


"CO" 
shell 


10^^ rad/s 


10-2 


10^^ erg X s 


He 
shell 


10~'^ rad/s 


10-3 


10^^ erg X s 


H 
shell 


lO-'^ratZ/i 


lo-*" 


lO^^erg X s 



4-3. The shape of rotating star 

Calculations presented above are based on the assumption that the den- 
sity stratification is not far from spheroidal and the rotation is close to 
"shellular", i.e. the angular frequency is constant on a family of spheroidal 
surfaces. The set of stellar structure equations is modified due to rotation 
but this method preserves the conventional description of the stellar struc- 
ture in terms of a one-dimensional mass coordinate, with constant mass 
surfaces being non-spherical [2H- In this approach some processes, such 
as the transport of angular momentum, the mixing of stellar matter and 
nuclear reactions are well described. However, the real two-dimensional 
structure is treated only approximately. 



Such an approach is especially unsatisfactory when we attempt to ana- 
lyze the collapse of the core because the star model is not in full mechan- 
ical equilibrium. The situation is analogous to starting the simulations of 
the rotating core collapse from the initial configuration which is that of a 
non-rotating core endowed with some amount of rotation. This case has 
been analyzed in [2b\ . where it was shown that some of models without 
initial equilibrium evolve in a significantly different way, compared to ini- 
tially hydro-stationary models. In this way only rough estimates of the basic 
quantities, such as e.g. the density at the core bounce, can be obtained. The 
analysis has been performed in the axisymmetric case with the cylindrical 
rotation and the barotropic EOS. One may expect a similar result for an 
arbitrary rotation law of the initial model. 
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Fig. 2. An example of the iso-density contours for n=3/2 (p = Kp^l^ 
EOS) polytropic model with a cylindrical distribution of the angular frequency, 
n(r) = r2o/(l + concentrated near the rotation axis (A = 0.1 i?Oj where i?o 

is the equatorial radius) calculated using eq. |5J|. The off-center density maximum, 
the flattened shape, the toroidal density distribution in the central region and the 
"cusp" in the polar region are common features of rapidly and differentially rotating 
polytropes. 

An important conclusion from the above considerations is that the exact 
density distribution and the velocity field inside the presupernova star in the 
full hydro-stationary equilibrium are required as an input for the analysis 
of the collapse of the core and its later evolution. To find the exact density 
distribution, in particular the shape of the stellar surface and the iso-density 
contours inside the star, one has to solve the equations of hydrostatic, or 
more explicitly hydro-stationary, equilibrium, taking into account the non- 
zero velocity field. This has been achieved in very simplified stellar models 
only. For the equation of state p = p(/o), including polytropes, the rotation 
was assumed to be cylindrical with angular velocity being a function of the 
radial variable in cylindrical coordinates only j^. Under these conditions 
one can derive an analytic first-order approximation to the resulting integral 
equation for the axisymmetric density distribution: 



pi(r,z)=rM/(po)-'J>c(r) + ($c>], 



(2) 
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where p is the density as a function of cyhndrical coordinates (r, (j), z), and 
($c) is the centrifugal potential <I>c corresponding to the angular frequency 
distribution inside the star r2(r), 



Jo 

averaged over the volume Vq of a non-rotating star with the density po{r) 
and the radius Rq: 



In eq. © f{p) is the enthalpy^ of the barotropic gas with an EOS p = p{p): 



From eq. (or using any numerical method of Ref. j261 127j ) we can see 
(Fig. [21) that the shape of the stellar core is far from spherical or spheroidal 
one. In case of the rapid differential rotation the density distribution is 
toroidal with an off-center maximum. This "disk-like" shape (Fig. [21) is 
sometimes referred to as a "concave-hamburger". A common feature of 
realistic density calculations is a cusp at the star surface in the polar region. 
It is a result of the strong centrifugal force near the rotation axis which is 
balanced by a rather weak gravity far from the center. 

Using the equilibrium approach we can obtain highly deformed and very 
flattened structures, but we can not forget that if they exceed stability lim- 
its their structure will be destroyed. Less constrained rotation, different 
from cylindrical and rigid one, results in a much more complicated prob- 
lem because a physically consistent treatment requires accounting for the 
temperature-dependent equation of state and the meridional circulation in 
a star. The importance of inclusion of the "shellular rotation" and the ex- 
ternal pressure is difficult to asses at present. The conclusion is that the 
current status of the stellar structure models with rotation does not allow 
us to predict firmly the initial state of the collapsing core in the supernova 
theory. A non-axisymmetric and ring-like shape of the core is a very intrigu- 
ing possibility as it may be partially responsible for violent processes during 
the last phases of the presupernova evolution and for the onset of the core 
collapse. The growing evidence of both axisymmetric and non-axisymmetric 
supernova explosions indicates a significant amount of rotation of the presu- 
pernova core but details are still unclear. It seems thus quite possible that 
at least some of the supernova progenitor cores can be properly described 
within a framework of the sketched above picture. 

^ The density in eq.(2) is the inverse function f-^: .f{f~\0) = f~\.f{C)) = C- 




(3) 






(5) 
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4.4- Rotating core collapse 

We mentioned already that the ratio Erot / Egrav grows as a result of the 
core contraction. This happens during subsequent nuclear burning phases 
and collapse. When the star exhausts a given nuclear fuel (H, He, C/0, Si) 
the core contracts until the next nuclear reactions ignite. After contraction 
the size of the core is smaller but the ratio Erot/ Egrav is bigger. During last 
phases of nuclear burning it may reach values as high as 0.04 j24j . 

To imagine what may happen in the collapsing and rotating core let's 
look at the behavior of a homogeneous and rigidly rotating sphere of a 
decreasing radius. In this case: 

^rot -J _ ^rot fM) 
Egrav EfXv ^ 

where quantities with indices "0" are just before the beginning of shrinking. 
From eq. © we can see that the ratio Erot/ Egrav grows as ^R~^. From 
this simple consideration one can infer the existence of the critical radius 
for which every stability limit will be exceeded. This implies impossibility 
of too big a shrinking even for a moderately fast rotating body. The process 
of shrinking will be stopped by the body breakup or another catastrophe. 

During the collapse of the core its radius, R, is decreasing very fast. 
As we mentioned, in calculations of ^24^ the iron cores with Erot/ Egrav = 
0.01 . . . 0.04 and with the initial radii Rq ~ 2000km are found. The final ra- 
dius is « 10km - a typical value of the neutron star radius. Actually, just 
after the collapse the radius of a newborn protoneutron star is -R ~ 60km, 
and after about one second it shrinks to ~ 20km. Finally, after tens of 
seconds the radius sets at R ^ 10km as a result of deleptonization. The 
dynamical instability limit of 0.27 for Mac Laurin spheroid is exceeded if 
the Fe core radius shrinks to 80 . . . 300km, respectively, for the initial ratio 
Erot/ Egrav = 0.01 . . . 0.04. This indicates the possibility of a violent hydro- 
dynamical instability. Because the shock front is born at approximately the 
same radius and at the same time, only very detailed dynamical simulations 
can give the actual behavior of the collapsing core. 

Recently, results of full three-dimensional simulations with no symmetry 
of the collapsing core assumed have been published [25]. In this paper a 
simplified equation of state and initial models of US] were used. This work 
focused on the gravitational radiation emission from the core. The core gets 
rid of the extra angular momentum in the form of the spiral arms or a small 
"satellite" as can be seen (see footnote |7| on page El for WWW address) 
in the results of Rampp [^. Possible future detections of gravitational 
waves correlated with SN explosions by the next generation detectors (LISA, 
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LIGO) will be an ultimate proof of the strongly asymmetric processes in the 
supernova engine. 

Detailed numerical calculations of the supernova explosion with inclusion 
of the neutrino processes have been reported in '29' '301 for an axisymmetric 
case with the equatorial symmetry assumed. However, the obtained results 
are rather ambiguous (see next section). 

5. Rotation and extended classification scheme 

In previous sections we have discussed some problems encountered when 
one attempts to include rotation in the core-collapse supernova theory. A 
question thus arises is such an extension of the theory necessary from the ob- 
servational point of view? The answer is definitively affirmative. Presently 
available astronomical data relevant to the supernova phenomenon suggest 
strong non-sphericity of the explosion. The most important pieces of evi- 
dence include velocities of pulsars sometimes correlated with an apparent 
deformation of remnants and the asymmetry of explosion deduced from the 
measured polarization. Also, some events appear to be superluminous under 
assumption of the spherical symmetry. Finally, the very recent identifica- 
tion of the supernova component in the light curve of GRB011121 provides 
a compelling argument in favour of the core rotation as discussed below. 

Let us explain how the core rotation can account for the those features 
of the data that can not be understood in the standard supernova scenario. 

Pulsars have much higher velocities than the average velocity of stars in 
the disk of the Galaxy. The highest observed speed exceeds 2300 km/s 
The observed pulsar speed [32 is generally a sum of two components: 

• a "kick" imparted to the nascent neutron star during the supernova 
explosion 

• an inherited velocity due to progenitor revolution in a binary system 
- when SN occurs in a binary (what happens quite often, probably 
most SN are in binaries) 

If the explosion disrupts a binary system we can see a runaway pulsar. The 
existence of binary pulsars tells us that some of systems survive two SN 
events. This also generally requires some momentum transfer to reduce the 
orbital velocity because of dramatic mass loss during the explosion. The 
velocity component due to the pulsar's birth in a binary system can be 
subtracted statistically from the distribution of pulsar velocities. We then 
obtain the distribution of kick velocities directly related to the SN explosion 
(Fig. 121). These velocities are still large enough to challenge the results of 
spherically symmetric supernova models. Proposed explanations of such 
large kicks include |34j : 
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Fig. 3. Distributions of pulsar velocities. The solid line corresponds to the ob- 
served velocities and the dashed line to natal velocities obtained by subtracting 
binary effects. Both solid and dashed curves are described by a x^-like distribu- 
tion, f{vp) — avp^^ exp{—Vp/b), adopted from with a = 1.96 x 10^® b = 514/3 
and a — 2.7 x 10^^ b — 60, respectively. 

• hydro dynamically driven kicks (fluctuations in a non-rotating core col- 
lapse) 

• asymmetric neutrino emission 

• electromagnetic post-natal "rocket effect" 

All these models disregard any dynamically significant rotation. As a kine- 
matical effect, the rotation adds some stochastic component to such pro- 
cesses as the convective motion or the neutrino emission. 

Although kicks in rapidly rotating supernova models vi^hich lead to non- 
axisymmetric deformations have not been studied systematically, the asym- 
metric disruption of a rapidly rotating core can be seen in some calculations. 
Results presented in 2_8_ show how a non-zero component of the velocity 
perpendicular to the rotation axis arises in simulations of a formation of 
PNS^. 

It is in principle possible to distinguish between kicks produced with 
a different amount of rotation if the relative orientation of the spin and 
the direction of motion is known. When the rotation is slow with respect 

^ http: / /www. mpa-garching. mpg.de / ~ wfk/MO VIES / 
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to the kick mechanism there is no correlation between the rotation axis 
and the kick direction. When the rotation is dynamically unimportant 
but faster than the kick mechanism, the angular averaging effect leads to 
the parallel orientation. For very fast rotation, leading to the asymmetric 
disruption of the core, we expect the rotation axis to be perpendicular to 
the momentum. In all three cases the relative orientation of the current 
velocity vector and the rotation axis can be found. The analysis of the 
binary neutron star system B1913+16 fH3| gives the result that the angle 
between the kick direction and the orbital plane is less than 5° — 10°. If 
the orbital momentum and the spin were aligned before the explosion, as 
it was assumed in the cited article, we would get an almost perpendicular 
orientation of the kick and the PNS rotation axis. The study of Vela and 
Crab pulsars [HII based on observed jets shows an apparent alignment on 
the sky of the velocity and the rotation axis. 

Statistical analysis of other 28 pulsars shows no correlation between their 
motion and rotation • Ambiguity of these results is due to difficulties in 
the analysis itself, simplifying assumptions and poor quality of the data. We 
only note that perpendicularity of the rotation axis and the kick, favored by 
the asymmetrical disruption of a rapidly rotating core in the SN explosion, 
is not disproved from the observational point of view. 

We also note that the fast rotation of pulsars is a strong evidence in 
favour of a huge angular momentum of the central part of exploding star. 
Rotation of the neutron star is the immediate consequence of the presuper- 
nova core rotation. There is no need to invoke any other mechanism, such 
as e.g. an off-center kick, to account for fast rotation of pulsars, in contrast 
to stellar models with very slow rotation of the presupernova core. 

5.1. GRB - SN connection 

Next argument in favour of strong rotation in supernovae comes from 
observations of superluminous events like SN1998bw - well known from its 
controversial coincidence with the gamma ray-burst GRB980425. Spectral 
analysis classified it as type Ic, but the absolute magnitude was close to that 
of type la j36ll37j . A spherically symmetric modelling has given an enormous 
explosion energy of 20-50 foe, the ejected mass in the range 12 — ISM© and 
the radioactive nickel amount as big as 0.5 — O.SMq |38j . Listed values are 
considered impossible in the standard supernova scenario. That is why these 
events are often called hypernovae^. As it has been shown in [HH] observed 
results could be well understood if the prolate asymmetry of the expansion, 
with the axis ratio 2:1, is assumed. We note, that the prolate expansion 



This in not a hypernova in the sense of the theoretical GRB models, but the possibility 
of a real connection between the two is not excluded and is a topic of current research 
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velocity produces oblate iso-density contours. This requires the standard 
explosion energy of 2 foe if observed 60° above the plane of symmetry. 
The ejected mass and the ^^Ni amount are 2Mq and 0.2Mq, respectively 
- typical values for a luminous core-collapse supernova. Polarization of 
light, which is more directly related to the asymmetry of expansion, is often 
present during hypernova events, again indicating a non-spherical explosion. 
Analysis of this sort, based on the rescaling of a spherical model tells us 
nothing about physical processes responsible for the asymmetry, but the 
rotation of the collapsing core is one of the most probable reasons. 

Very recently, observations of a SN event in the light curve of GRB011121, 
have been reported 10 . Joint campaign of optical and X-ray observations, 
also from space, succeeded to firmly establish the presence of a supernova 
in the GRB011121 place, when the afterglow has declined. The supernova 
is of type Ic. Collected observations of GRB011121 fit nicely to the collap- 
sar model of MacFadyen and Woosley |41j. In this model, rotation of the 
progenitor's core plays a crucial role and the rotation axis determines the 
direction of jets which produce observed gamma rays. Despite the formation 
of a black hole from the collapsing core, matter rich in ^^Ni is expelled along 
the rotation axis and its emission forms a typical supernova light curve ||4lj| . 

5.2. Including rotation into the classification scheme 

We feel encouraged by the above analysis to unify the phenomena of su- 
pernovae, hypernovae and some GRB's in a two-dimensional classification 
scheme which will expand that shown in Fig. ^ As a second dimension 
we include some measure of rotation of the progenitor's core {e.g. its an- 
gular momentum). We conjecture that classical supernovae correspond to 
some relatively low amount of angular momentum, which is just enough 
to produce a neutron star with proper natal momentum. Here we assume 
that the explosion itself may be treated as a strong argument in favour of 
a sufficiently high rotation of the presupernova core. The simulations have 
shown the failure of the idea of explosion by the prompt mechanism. The 
delayed mechanism improved the situation but it is not clear at all if it 
really works. Indeed, some of the models still fail to produce SN events in 
spite of inclusion of the neutrino processes. 

When the core angular momentum if bigger one expects that some part 
of the collapsing core is locked for a while in a sort of "fizzler" and is 
accreted back onto the inner part. Depending on the falling back mass, a 
neutron star can form, or further collapse to black hole can proceed. As 
the fizzler possess large angular momentum, very fast rotating neutron star 
can form, if the accreted mass does not exceed some fraction of the solar 
mass. In case the mass is much higher, of a couple of solar masses, as 
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Fig. 4. Two-dimensional classification of supernovae and some GRB events. The 
vertical axis represents a property which is can serve as a measure of the core 
angular momentum 



in numerical simulations happens for heavier iron cores formed in massive 
helium stars (of say 10 solar masses), a Kerr black hole finally forms. A 
transient accretion disk feeds angular momentum to the black hole and 
cools emitting strong neutrino flux. One can conjecture, that hypernovae 
and GRB's are displays corresponding to such a way of the core collapse. 
Positions of various phenomena in the two-dimensional classification scheme 
are shown in Fig. ^ 



5.3. Discussion 

If simulations of the rotating supernova fail in the sense that the energy 
of radiation and ejecta is similar to that obtained in the simulations without 
rotation which give a sufficient explosion power we have to accept standard 
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scenario. In the opposite situation we have an indirect argument in favour 
of rotation of the Fe core. The reasoning presented above holds provided the 
results of simulations of such a complex problem as the supernova engine are 
reliable. One should notice that some serious difficulties are reflected in the 
published results of simulations with rotation. Yamada and Sato [HO] found 
stronger explosion but recently Fryer and Heger j29j obtained decrease of 
explosion energy. 

Finally we note that the evidence of strong asphericity in supernova 
explosions is being established now. Besides the arguments mentioned in 
the previous section, some very suggestive images of supernova remnants 
are available. One of them is Puppis A with the asymmetry apparently 
correlated with the movement of a young neutron star |42|. Pictures like 
this are very suggestive and tell us that pulsar kicks are natal and produced 
by the supernova engine. A lot of processes may be responsible for this, 
however rotation is the most common feature of astrophysical objects and it 
is reasonable to examine its possible effects first. The knowledge relevant to 
a wide range of problems related to rotation, such as the initial pre-collapse 
state and the behavior of the contracting core which exceeds the rotation 
stability limits, is still unsatisfactory. Advances in the supernova rotation 
research will hopefully lead to a better understanding of many problems 
of supernova physics, such as the mechanism of pulsar kicks, the GRB-SN 
connection, the emission of gravitational waves and the formation of Kerr 
black holes. 

The authors are grateful to W. J. Swiatecki for critical reading of the 
manuscript and for claryfing discussion on stability of Jacobi ellipsoids j43j . 
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